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ABSTRACT 

Cirrus cloud radiative and physical characteristics are determined using a combination of ground-based, 
aircraft, and satellite measurements taken as part of the FIRE Cirrus Intensive Field Observations (IFO) during 
October and November 1 986. Lidar backscatter data are used with rawinsonde data to define cloud base, center, 
and top heights and the corresponding temperatures. Coincident GOES 4-km visible (0.65 pm) and 8-km 
infrared window ( 1 1.5 (im) radiances are analyzed to determine cloud emittances and reflectances. Infrared 
optical depth is computed from the emittance results. Visible optical depth is derived from reflectance using a 
theoretical ice crystal scattering model and an empirical bidirectional reflectance model. No clouds with visible 
optical depths greater than 5 or infrared optical depths less than 0.1 were used in the analysis. 

Average cloud thickness ranged from 0.5 km to 8.0 km for the 71 scenes. Mean vertical beam emittances 
derived from cloud-center temperatures were 0.62 for all scenes compared to 0.33 for the case study (27-28 
October) reflecting the thinner clouds observed for the latter scenes. Relationships between cloud emittance, 
extinction coefficients, and temperature for the case study are very similar to those derived from earlier surface- 
based studies. The thicker clouds seen during the other IFO days yield different results. Emittances derived 
using cloud-top temperature were ratioed to those determined from cloud-center temperature. A nearly linear 
relationship between these ratios and cloud-center temperature holds promise for determining actual cloud-top 
temperatures and cloud thicknesses from visible and infrared radiance pairs. 

The mean ratio of the visible scattering optical depth to the infrared absorption optical depth was 2.13 for 
these data. This scattering efficiency ratio shows a significant dependence on cloud temperature. Values of mean 
scattering efficiency as high as 2.6 suggest the presence of small ice particles at temperatures below 230 K. The 
parameterization of visible reflectance in terms of cloud optical depth and clear-sky reflectance shows promise 
as a simplified method for interpreting visible satellite data reflected from cirrus clouds. Large uncertainties in 
the optical parameters due to cloud reflectance anisotropy and shading were found by analyzing data for various 
solar zenith angles and for simultaneous AVHRR data. Inhomogeneities in the cloud fields result in uneven 
cloud shading that apparently causes the occurrence of anomalously dark, cloudy pixels in the GOES data. 

These shading effects complicate the interpretation of the satellite data. The results highlight the need for 
additional study of cirrus cloud scattering processes and remote sensing techniques. 


1. Introduction 

Accurate quantification of cirrus cloud properties 
from satellite measurements is particularly important 
to the understanding of the role of cirrus in climate 
change. The nonblackness of cirrus at thermal infrared 
wavelengths renders the interpretation of satellite data 
taken over cirrus more difficult than measurements 
over most water clouds. The International Satellite 
Cloud Climatology Project (ISCCP; see Schiffer and 
Rossow 1983) is making an ambitious effort to derive 
daytime cirrus coverage, altitudes, and optical depths 
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over the globe during a 5-year period. The ISCCP anal- 
ysis algorithm (Rossow et al. 1988) relies entirely on 
bispectral data taken at visible (VIS: ~0.65 ^m) and 
infrared (IR: — 1 1.5 p.m) wavelengths. Although VIS- 
IR bispectral techniques have been suggested as feasible 
methods for determining bulk cirrus properties (e.g., 
Shenk and Curran 1973; Reynolds and Vonder Haar 
1977), there has been very little application of these 
techniques to real data prior to the ISCCP. 

The basic premise for using the bispectral approach 
is that the VIS extinction coefficient is related to the 
IR absorption coefficient. This relationship implies that 
the cloud VIS reflectance may be used to infer the 
cloud’s IR emittance. Having a value for the clear-sky 
IR radiance, it is possible to correct the observed cloudy 
radiance for cloud emittance resulting in an estimate 
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of the radiance emanating from a specified level in the 
cloud. The equivalent blackbody temperature of this 
level, usually the cloud center, is then converted to 
cloud altitude using a vertical sounding of temperature. 
The critical relationship required for this approach is 
the dependence of IR emittance on VIS reflectance 
through the IR and VIS optical depths. Since clouds 
scatter radiation anisotropically, this relationship is also 
influenced by the viewing and illumination conditions. 

The ISCCP cirrus analysis ( Rossow et al. 1988 ) uti- 
lizes a combination of theoretical and empirical models 
to determine the cloud visible optical depth from the 
observed reflectance; the cloud emittance from the vis- 
ible optical depth; and finally, the cloud-top temper- 
ature from the cloud emittance and the observed in- 
frared radiance. The theoretical cloud model is a ra- 
diative transfer scheme that simulates the scattering 
and absorption of visible radiation by water droplets 
with an effective radius of 10 ura. For water droplets 
of this size, the ratio of VIS extinction to infrared ab- 
sorption optical depths is ~2.4. An analysis of coin- 
cident satellite and lidar data by Platt et al. ( 1980) and 
theoretical calculations employing cylinders (Platt 
1979) suggest that this ratio is approximately equal to 
2.0 for cirrus. The ISCCP algorithm utilizes the latter 
value to provide a link between the water droplet model 
and actual cirrus clouds. 

Cirrus clouds are primarily composed of ice crystals 
with various shapes having maximum dimensions 
ranging from about 20 jxm to 2000 (e-g-> Heyms- 

field and Platt 1984). The scattering properties of hex- 
agonal ice crystals differ considerably from spherical 
particles (Liou 1986). Because of the complexities in- 
volved in computing scattering by hexagonal solids, 
cylindrical columns have been used to approximate 
hexagonal crystals in radiative transfer calculations 
(e.g., Liou 1973). More recently, however, Takano and 
Liou ( 1989a) have solved the radiative transfer equa- 
tions for randomly oriented hexagonal plates and col- 
umns. Their results are the most realistic to date in 
that they reproduce certain well-known cirrus optical 
phenomena. 

Absorption plays the dominant role in IR extinction 
in cirrus clouds. Some theoretical investigations ( Liou 
and Wittman 1979; Stephens 1980), however, have 
shown that scattering effects may also be significant at 
IR optical depths greater than ~0.1. An IR radiance 
measured by a satellite over cirrus clouds, therefore, is 
the product of both absorption and scattering processes 
in the cloud, as well as the transmission of radiation 
from below the cloud (Platt and Stephens 1980). It is 
generally assumed, however, that scattering effects are 
negligible so that the observed emittance is considered 
to be the absorption beam emittance. 

Empirical studies have also shed some light on the 
VIS reflectance-IR emittance relationship. Platt ( 1973 ) 
developed techniques for deriving cloud visible and 
infrared properties from a ground-based lidar and an 


upward-looking infrared radiometer. The backscattere 
intensities measured with the lidar are used to defin 
cloud base and top heights. Cloud emittance was dc 
rived from the observed down welling IR radiance. Plai 
and Dilley (1979) presented emittance results from 
set of observations taken over Australia. Platt et a 
(1980) used lidar and satellite VIS-IR data to estimat 
the dependence of beam emittance on VIS cloud re 
flectance for a limited set of viewing and illumination 
conditions over Colorado. Their results are more cor 
sistent with theoretical scattering from ice cylinder 
than with scattering from ice spheres. Aircraft radic 
metric measurements taken over New Mexico (Pa) 
tridge and Watt 1981) have also been used to determin 
the radiative characteristics of cirrus clouds as relate- 
to the cloud ice water path. Those results provide fui 
ther evidence that real clouds scatter more like cylinder 
than spheres. Platt (1983) combined the results fror 
previous studies and used them to explain the char 
acteri sties of two-dimensional bispectral histograms c 
VIS-IR data observed from a geostationary satellite 
Theoretical calculations of reflectance and emittanc 
for typical cirrus clouds were consistent with the sal 
ellite data taken over areas of suspected cirrus cloud: 
While that study provided encouragement for using 
bispectral approach to retrieving cirrus properties fror 
bispectral data, it also highlighted some of the diffi 
culties that are likely to be encountered with such 
technique. Platt and Dilley ( 1984) used lidar and sola 
radiation measurements to measure part of the single 
scattering phase function of real cirrus clouds. Thei 
results fell within the range of laboratory measurement 
and theoretical calculations for hexagonal crystals. A: 
analysis of a large sample of ground-based lidar an- 
infrared data taken over Australia (Platt et al. 1987 
showed that the average emissivity of cirrus clouds i 
primarily a function of the midcloud temperature 
Though fraught with significant uncertainties, tha 
study also indicated that the theoretical value of th 
ratio of visible extinction to infrared absorption fc 
cirrus clouds may be too low. 

From these previous studies, it appears that: 

1 ) cirrus cloud scattering properties are similar t 
those of hexagonal crystals resulting in reflectance pal 
terns that are unlike those from spheres; 

2) scattering of IR radiation may be important i 
determinations of IR optical depths; and 

3) the ratio of VIS extinction to IR absorption coel 
ficients is between — 1 .8 and 4.0. 

The full impact of these results on using a VIS-IR bi 
spectral mc,l.od for retrieving cirrus properties is un 
known. Differences between ice crystal and wate 
droplet bidirectional reflectance patterns will introdur 
errors into the retrieved VIS optical depth. Use of a 
observed beam emittance with a theoretical model tha 
assumes absorption only may affect the emittance e: 
timation. Finally, uncertainties in the extinction rati 
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(scattering efficiency) may cause significant errors in 
the estimation of IR optical depth. 

In this paper, the relationship between VIS reflec- 
tance and IR emittance is examined using data taken 
during the First ISCCP Regional Experiment (FIRE) 
Cirrus Intensive Field Observations (IFO; see Stan- 
1987). Ground-based and aircraft lidars are used to 
define the vertical locations of the cirrus clouds, while 
satellites provide measurements of VIS and IR radi- 
ances. Both VIS and IR optical depths are computed 
from the reflectance and emittance data covering a 
range of solar zenith angles missed in previous studies. 
These relationships are derived to provide a means for 
the application of a bispectral cirrus parameter retrieval 
algorithm over the FIRE IFO region. Results are pre- 
sented for the entire IFO period with emphasis on 27- 
28 October 1986, the case study period. The data pre- 
sented here also constitute an initial source for devel- 
oping cirrus bidirectional reflectance models and may 
be used to help validate the models employed in the 
ISCCP algorithm. 

2. Data 

a. Lidar measurements 

Lidar backscatter data were taken from four different 
sources — three surface and one airborne. The lidars 
and their operating systems and data products have 
been described elsewhere. Thus, only a brief description 
of sources and their uses in this study are given here. 

The University of Utah mobile polarization lidar 
(see Sassen et al. 1990) was located at Wausau, Wis- 
consin (WAU; 45.0°N, 89.7°W). The NASA Langley 
ground lidar (Sassen et al. 1990) was situated at Ft. 
McCoy, Wisconsin (FMC; 43.9°N, 90.8°W), while the 
University of Wisconsin High-Spectral Resolution Li- 
dar (Grund and Eloranta 1990) was in Madison, Wis- 
consin (MAD; 43.1°N, 89.4°W). These ground sys- 
tems acquired nearly continuous lidar backscatter pro- 
files during the cirrus days of the IFO with especially 
good coverage during the case study period. The lidar 
returns are used to define the cloud base and physical 
thickness. Under conditions of small attenuation and 
constant backscatter-phase function, the backscatter 
intensity profiles indicate the vertical distribution of 
cloud extinction. A time series of these lidar returns 
shown in Fig. 1 define the outlines of the cirrus clouds 
as they passed over FMC during the afternoon of 28 
October. Solid black areas define the most intense cloud 
backscatter. Grey denotes less backscatter and white 
indicates no cloud. Vertical stripes represent missing 
data. Cloud-top altitude is fairly constant at • — -11 km. 
Cloud base changes from ~8 km to 10 km approxi- 
mately every half hour. In the morning, cloud base was 
observed at ~7 km, while cloud top varied between 8 
km and 1 1 km. A similar variation is also seen in the 
WAU estimated volume backscatter coefficients shown 
in Fig. 2 for the morning of 28 October. Darker portions 



Fig. 1. Lidar backscatter ratios from cirrus clouds passing over 
Ft. McCoy during 28 October 1986. 


of the plot correspond to higher backscatter coefficient 
values. The particle backscattering efficiency depends 
on cloud particle shape and phase. Further details of 
the lidar returns are reported in the cited references. 

Three parameters are derived from plots like those 
in Figs. 1 and 2 by averaging the data within ± 1 5 min 
of the UTC (Coordinated Universal Time) half hour 
plus 5 min. All times, however, will be given here to 
the nearest half hour. Goud-top altitude, z t , and cloud- 
base altitude, z b , are defined as the average altitudes 
of the highest and lowest nonclear-air backscatter re- 
turns, respectively. Similarly, the cloud thickness is h 
= z, — z b . Mean cloud height (approximately cloud 
center height), z c , is the backscatter-intensity weighted 
average height of the cloud. It corresponds roughly to 
the altitude below which 50% of the lidar backscatter 
is accumulated. These parameters were estimated 
graphically for the FMC and WAU sites, while a com- 
puter analysis was applied to the MAD results. The 
value of Zc for MAD corresponds to the midpoint in 
optical thickness independent of cloud attenuation (see 
Grand and Eloranta 1990). Since the clouds are ad- 
verting over the fixed surface sites, the averaged lidar 
data correspond to a thin vertical cross section taken 
out of some cloud volume. It is assumed that the cross- 
section-averaged data represent the mean conditions 
of that volume. 

These same parameters were also derived from the 
down-looking lidar backscatter plots reported by Spin- 
hime et al. ( 1 988 ) for selected flight tracks of the high- 
flying, NASA ER-2 aircraft over the IFO area. Shorter 
time averages were used since the plane’s motion 
greatly increased the cirrus advection rates relative to 
the lidar. In some instances, the clouds were too thick 
for complete penetration by the ground-based lidars. 
To determine these occurrences, the cloud altitudes 
estimated from the ground were compared to those 
determined from the nearest aircraft flight. On most 
days, there was good agreement between the surface 
and airborne lidars. The thick clouds observed on 22 
October required use of the aircraft lidar to estimate 
z,. At other times when no aircraft data were available 
for comparison, a different approach was used to es- 
timate z, (see section 3a). 
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Fig. 2. Estimated lidar volume extinction efficiencies 
from cimis clouds over WAU 28 October 1986. 


b. Satellite radiances 

Half-hourly VIS and IR data from the Geostationary 
Operational Environmental Satellite (GOES) were re- 
solved to 4-km pixels by averaging and replication of 
1-km VIS and 4-km X 8-km IR pixels, respectively, 
for most times. When the full-resolution data were not 
available, 8-km pixels derived by pixel averaging ( VIS ) 
and sampling (IR) were used. The pixels were navi- 
gated to ±8 km as in Minnis and Harrison ( 1984a). 
Pixel pairs were extracted for areas within 0.35° lon- 
gitude and 0.25° latitude of the surface lidar sites. These 
areas correspond to cloud advection at wind speeds of 
~ 30 m s' 1 for a half-hour interval. Strips of pixels, 
4(2) pixels wide, along the wind vector at z c centered 
over the surface site were then taken from the 4 km( 8 
km) resolution areas. The satellite data corresponding 
to the ER-2 flight track were extracted for a strip cen- 
tered on the flight track. This strip is approximately 
30 km wide and 100 km long. 

Two-dimensional histograms of the VIS and IR ra- 
diance pixel pairs were formed from the strips of ex- 
tracted data. The VIS data, stored as counts, D (where 
0 < D « 63), in the histogram were converted to ra- 
diance and VIS reflectance, p, using the calibration of 
C. H. Whitlock ( 1989, personal communication) where 

P = (0.1624Z) 2 - 8.3)(526.9m 0 )-*, 

and po is the cosine of the solar zenith angle $ 0 . The 
raw IR data are given as equivalent blackbody tem- 
peratures, T, and are converted to radiance with the 
Planck function, B(T), evaluated at 11.5 pm. Note 
that all radiometric quantities discussed in this study 
are spectral quantities, either VIS or IR, so no subscripts 
denoting spectral dependence are given. 

A similar set of histograms was derived from the 


NOAA-9 Advanced Very High Resolution Radiome 
(AVHRR) Global Area Coverage (GAC) 4-km d 
for the 0.35° X 0.25° areas over the surface sites. 1 
NOAA-9 daytime pass over the areas occurred betwe 
2000 and 2130 UTC ( — 1430 local time). Radiam 
from channels 4 and 5 were averaged to produce 
IR radiance corresponding to the GOES 11.5-p 
channel. The AVHRR VIS data were taken fre 
channel 1 and converted to reflectance using the ca 
brations of C. H. Whitlock ( 1989, personal comm 
mcation) where 

P = (0.6060 D a - 22.0)(519.4po) -1 

and D a is AVHRR 10-bit count. 

Solar zenith, satellite zenith 6, and relative azimu 
Jp an £lcs were computed for each set of measuremen 
from GOES-6 located over the equator at appro) 
mately 97°W. This location yielded a value of 8 ~ 
for the ground sites. Relative azimuth angles were coi 
fined to the backscattering hemisphere for GOES. T1 
orbital precession of the NOAA-9 produces a variatic 
in 0 between 0° and 70° over a given site every 4- 
days. Because of its cross-track scan pattern, th 
AVHRR views a given site from relative azimuth anglt 
that lie in both the forward and backward hemisphei 
at a nearly constant skew to the solar plane. The ske' 
depends on both the solar and site latitudes. Const 
quently, a range of solar zenith angles with constar 
viewing anHes is covered using the GOES, while a rang 

of viewing angles at a relatively constant do is covere 
with the AVHRR. 

Days when cirrus were observed over the IFO with 
out substantial low-cloud interference are referred t< 
as cirrus days. The times used in this study when lida 
and satellite data coincided with cirrus are listed ii 
Table 1. 
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Table I. Times and locations of lidar-satellite data 
used in this study. 


Site 

Day 

Month 

Times (UTC) 

FMC 

22 

October 

1300, 1330, 1400, 1600, 1630, 1700, 




2000 


27 


2030, 2100 


28 


1330, 1400, 1430, 1500, 1600, 1700, 
1900 

1930, 2000, 2030, 2100, 2130, 2200 


30 


2000, 2030 


1 

November 

1800, 1900 


2 


1900,2000,2100 

MAD 

28 

October 

1330, 1500, 1600, 1700, 1800, 1930, 
2000 

2030,2100,2130, 2200 

WAU 

22 

October 

1300, 1330, 1400, 1430, 1600, 1630, 
1700, 1800 

1830, 1900, 1930,2130, 2200 


28 


1500, 1600, 1700, 1800, 1900, 1930, 
2000, 2030 
2100, 2130 


30 


2130, 2200 


1 

November 

1800, 1900, 2000 


2 


1700, 1800, 1900,2000,2100 


c. Temperature data 

Soundings from Green Bay, Wisconsin, determined 
the temperature-height relationships for all of the data. 
Linear interpolation was used to estimate half-hourly 
soundings from the six-hourly data. Cloud-top tem- 
perature, T t , corresponds to z t on the soundings. Mean 
cloud temperature, T c , is found from z c . Surface tem- 
peratures taken every six hours at MAD, WAU, and 
Lone Rock, Wisconsin (Hahn et al. 1988), and oc- 
casionally at the FMC site, were used to supplement 
the clear-sky temperatures derived from the satellite 
data as described below. 

The clear-sky temperature, T s , is the equivalent 
blackbody temperature for clear scenes. It is estimated 
in several different ways. The first order estimate is 
taken from the initial results of Minnis et al. ( 1990), 
which applies the techniques of Minnis et al. ( 1987) 
to 0.5° regions within the greater IFO area. That ap- 
proach sets a VIS threshold — 2 counts above the clear- 
sky count, D t (see section 3). All pixels considered to 
be clear must be darker than this threshold and have 
a temperature that is no more than 3 K colder than 
the maximum observed temperature. The 4 K range 
for clear pixels, roughly double the typical value over 
Wisconsin land areas, allows for shading effects in 
partly cloudy scenes. The average temperature of the 
clear pixels is the initial value of T t . Surface air tem- 
peratures, T t , are also taken from nearby ground sta- 
tions. A rough correction is applied to these tempera- 
tures to adjust for atmospheric attenuation and the 
difference between the temperature of the surface skin 
and the air at shelter height. The resulting estimate of 


clear-sky temperature is T t . An example of the rela- 
tionship between T t and T t is shown in Fig. 3 for a 
region including FMC during 27 October 1986. Note 
that the difference between T, and T % constitutes the 
correction to determine T t from T t . It varies with local 
time (6 hours less than UTC) from positive values in 
the morning to negative values in the afternoon, be- 
havior typical of land surfaces (e.g., Minnis and Har- 
rison 1984a). 

The initial value of T s for the ground site is compared 
to the values for the surrounding regions to check for 
cirrus contamination of the clear-sky temperature. If 
r, for the ground site is not within 2 K of the warmest 
value found for nearby regions, then its value is reset 
to that of the warmest value. The 2 K limit is typical 
of the range in T t over relatively flat homogeneous 
areas. If extensive overcast prevents the calculation of 
T s , T s — T t . 

The AVHRR clear-sky temperatures are first esti- 
mated using T % derived from the AVHRR data. If that 
result is lower than the corresponding GOES estimate 
of by 2 K, then the AVHRR value of T % is reset to 
the GOES value. It is assumed that the warmer of the 
two values is less cloud contaminated. Because of limb- 
darkening effects, the AVHRR value of I , may be up 
to 2 K colder than the GOES value and still be valid. 
In warmer, moister atmospheres, greater limb-dark- 
ening effects are expected. 

3. IR and VIS parameter analyses and results 

a. Emittance calculations 

Neglecting IR scattering effects, the observed cloud 
beam emittance is given here as 



UTC 

Fig. 3. Comparison of clear-sky and shelter air temperatures over 
FMC 27 October 1986. 
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«b< 0 ) = [ 5(7 ) - 5(7,)][5(7 l ) - 5(7,)]-', (1) 

where 7 Z is the temperature at some altitude z corre- 
sponding to the cloud. The mean clear-sky equivalent 
blackbody temperature over the area of interest, 7,, 
has a weak dependence on 0. 

Cloud beam emittance is calculated twice for each 
set of lidar-radiance data using 7 Z = 7 C and 7 Z = 7,. 
The former value, which corresponds to the quantity 
used in most previous studies (e.g., Platt et al. 1980), 
may be more representative of the actual radiating part 
of the cloud. It does not necessarily correspond to the 
center of the cloud. The actual cloud-top temperature 
defines the vertical limit of the cloud. Emittances are 
computed using both temperatures to determine if 
there is a relationship between them that may be used 
to better define the physical boundaries of the cloud 
from the satellite data. 

Due to angular effects, the values of 7 C may require 
some adjustment from the initial lidar values. It is un- 
likely that a value of « b = 1 will be measured at a useful 
satellite zenith angle using 7, because of the low density 
of particles in the upper portions of the cloud. On the 
other hand, « b ( 7 C ) may be greater than one for some 
thick cirrus clouds. Although emittances greater than 
unity may be possible due to scattering enhancements 
of the upward radiance (Platt and Stephens 1980), the 
uncertainties in z c for thick clouds preclude any defin- 
itive measurements of ^ > 1. Thus, if e b > 1, 7 C is 
decreased until the average value of t b for a given re- 
flectance is less than or equal to one. While this limit 
is reasonable, it is somewhat arbitrary resulting in in- 
creased uncertainty in the true value of 7 C for thick 
clouds. 

During initial processing of the data, it was deter- 
mined that the maximum emittance found using 7, 
was ~0.86, except for those cases with cloud cover too 
dense for complete penetration of the lidar beam. To 
identify and correct the exceptions, a new estimate of 
T x was computed whenever 7 C was adjusted as ex- 
plained above. This new estimate, T[, is determined 
from the following formula, which forces the cloud to 
have the maximum observed emittance: 

5(7' t ) = (5(7) -0.145(7 z )]/0.86. 

The resulting value was then compared to 7 t . If T\ 
<T t - 3 K, then 7, is reset to T[. The 3 K allowance 
( — 0.5 km) is made to account for uncertainties in the 
lidar-determined cloud-top, due to time averaging and 
unknown penetration depth in thicker cloud. The result 
was then compared to the tropopause temperature us- 
ing the assumptions that the cloud occurs in the tro- 
posphere and the tropopause temperature is the coldest 
in the troposphere. If 7, is colder, it is reset to equal 
the tropopause temperature. The value of z, was then 
adjusted to correspond to the final value of 7 t . 

It is also assumed that 

e b = 1 - exp(-r e /p), (2) 


where r e is the IR absorption optical depth ar 
= cos0. Based on the results of Platt and Step! 
(1980), it is expected that the viewing zenith a 
dependence of e b will not depart significantly from 
Values of beam emittance derived with AVHRR 
may be adjusted to the GOES viewing zenith a 
with this relationship. 

The vertical emittance from (2) is 

«. = 1 - exp(-T e ). 

It is assumed here that scattering effects are neglij 
in the upwelling direction. Thus, t c is equivaler 
the IR absorption optical depth and t, is equal tc 
vertical emittance. 

b. VIS reflectance and optical depth calculations 

Values of clear-sky reflectance p, and clear-sky o 
D, were computed for each region using the 0.0 1 ° c 
sky albedo, a x , map of the IFO area (42°N-4. 
87°W-92°W) constructed by Minnis et al. ( 1990) f 
GOES data at each half hour. Clear-sky reflectance 
any latitude X and longitude 0 of the grid at time 
estimated as 

p s (X, 0, t, 0 O , 9, 0) = a t (X, 0, t, 0 O )X S (0 O , 0, 0), 

where X, is the anisotropic reflectance factor with v£ 
given by the model of Minnis and Harrison ( 198 
D t was determined by p, using the VIS calibrati 
The value of 0o varies by a few degrees over the 
time period, while the values of a t were normalize 
a single value of solar zenith angle designated 0o, 
account for these variations, a,(t, 0o) = < 
0 o )po/mo i > where po, = cos0o,. The clear-sky diffus 
bedo is 


«sd 


= J a,(0 o )poduo I J uoduo. 


integrated over po = 0, 1. The value of a ^ is set e 
to a,(57°) in this study since the full range of s 
zenith angles is not observed at the time and lati 
of the IFO, the true value of is usually equivj 
to a* measured at 0 O «= 53°, and 57° is the lo 
observed 0o for this dataset. 

Cloud reflectance, p c , is estimated with a variai 
the simple physical model used by Platt et al. ( 19 
That is. 


p - T,p c + p,T c T u + a*( 1 - Od)( 1 - 7 C - a*), 

where p is the measured reflectance, a c is the c 
albedo at 0 O , Xc is the anisotropic reflectance facto 
the cloud, p c = OcX c (0o, 0, 0). 

This model assumes that all ozone absorption oc 
above the cloud (first term ) and all Rayleigh and s 
sol scattering is confined to the layers below the cl 
The second term in (5) accounts for direct sola 
diation, which passes through the cloud, reflects f 
the surface, and passes back through the cloud ir 



November 1990 


MINNIS, YOUNG, SASSEN, ALVAREZ AND GRUND 


2408 


direction of the satellite. The third term accounts for 
the radiation that passes down through the cloud via 
multiple scattering, reflects diffusely from the surface 
below the cloud, and returns through the cloud scat- 
tered in the direction of the satellite. 

Using the parameterization of Rossow et al. ( 1988 ), 
the transmittance of the air above the cloud is 

T. = exp[-u(0.085 - 0.00052w)( 1/mo + 1/m)]. 

where u is the ozone abundance in cm-STP. The value 
used here, u = 0.32 cm-STP, is the average of the mid- 
latitude winter and summer standard atmospheres 
above 10 km from McClatchey et al. ( 1973). Platt et 
al. ( 1980) implicitly assumed that T a = 1 . The current 
model accounts for ozone absorption in the Chappius 
bands. 

The transmittance of the cloud to direct solar radia- 
tion at 6 0 is 

T c = exp(-r v /2/io), (6) 

(see Platt et al. 1980). Similarly, the direct transmit- 
tance from the surface through the cloud along the 
satellite line of sight is 

T u = exp(-T v /2M). 

The visible optical depth is reduced by a factor of two 
for the direct transmittance because at least half of the 
radiation scattered out of the beam is actually diffracted 
in the forward direction (Takano and Liou 1989a). 
Clear-sky reflectance along the satellite line of sight is 
p s and a,,) is the effective clear-sky albedo to diffuse 
radiation directly below the cloud. Due to the relative 
homogeneity of clear-sky reflectance over the IFO re- 
gion, it is assumed that and p s may be computed 
from the same data. The albedo of the cloud to diffuse 
radiation is a d . 

In addition to values for the clear-sky terms, the so- 
lution of (5) for p c requires specification of r v and X c . 
VIS optical depth is estimated by iteration on ( 5 ) using 
a linear interpolation of the relationships between po 
and a c for randomly oriented hexagonal columns 
(length, 125 pm; width, 50 pm) in Fig. 4 of Takano 
and Liou ( 1989b). Similar interpolations are used to 
estimate a d (r v ), where 

«d(r v ) = 

Jo 

For a given measurement, (5) is solved iteratively using 
an initial guess of cloud albedo such that T c = T u = 1 
— a c . A value for r v is determined from this initial 
guess using the theoretical data. A limit of 20 iterations 
is imposed to achieve an absolute difference of less than 
0.001 between the guess and the computed value of 
a c . Generally, fewer than five iterations are required. 
Since a c must be greater than zero, a c is set to 0.001 
for initial guesses that are less than or equal to zero. If 
t v < 0, a d < 0.00 1 , or a c 0.00 1 after any iteration. 


Od(r v , po ) dfio / Po^Po- 


it is assumed that r v is indeterminate and the data are 
not used. The causes for these indeterminate cases (e.g., 
shadowing of the surface by adjacent clouds or inad- 
equate specification of the cloud reflectance anisotropy ) 
are discussed later. 

A value for X c , which depends on r v and the cloud 
microphysics, is needed to determine a c from p c . No 
models of Xc are currently available for ice clouds in 
terms of r v . Because of favorable angles Platt et al. 
(1980) were able to assume that X c = 1. However, most 
empirical and theoretical bidirectional reflectance 
models for cloudy scenes (e.g., Suttles et al. 1 988 ) reveal 
a systematic decrease in x c with 0 O for the angles used 
in this study. The cloudy scene bidirectional reflectance 
model developed by Minnis and Harrison ( 1984b) is 
used initially to estimate Xc. That model’s reflectance 
anisotropy is similar to other empirical and theoretical 
models (Stuhlmann et al. 1985). The inclusion of all 
cloud types in its derivation should produce a reflec- 
tance pattern that combines the scattering properties 
of both ice and liquid water clouds. New values of X c , 
derived after the initial analysis, are used to reanalyze 
the data to provide better estimates of VIS optical 
depth. 

The data were preprocessed to define limits to elim- 
inate pixels containing low clouds and those scenes 
containing only thick clouds. Underlying low clouds 
confuse the interpretation of cirrus radiances. Thick 
clouds increase the uncertainty in the determination 
of T c , T t , and t b . A simple filter of the form, 

f b = 1 - exp(-fca c /p), (7) 

where k is a regression coefficient, was used to eliminate 
low clouds. This formula gives a first approximation 
to the relationship between t b and cloud albedo. The 
rationale for its use and the details of the filtering are 
described in appendix A. The data were also screened 
for partially cloud-filled pixels as detailed in appen- 
dix B. 

c. Results and discussion for midcloud temperature 
emittances 

All results discussed in this section are based on T z 
= T c in ( 1 ). Examples of the two-dimensional GOES 
histograms used in this analysis are shown in Figs. 4a 
and 4b for 1500 UTC over FMC. The latter represents 
a cirrus case (see Fig. 1 ) on 28 October, while the 
former, taken during the previous day, is typical of 
clear conditions. Maximum clear-sky reflectance for 
this hour is denoted with the dashed line in Fig. 4a. 
Some of the cold, apparently cloudy pixels in Fig. 4b 
are no brighter than the clear pixels in Fig. 4a. More- 
over, some of these pixels are actually darker than the 
cloud-free pixels. Depending on X c and r v , some of the 
cold, dim pixels yield a positive value of a c in the so- 
lution of (5). Those pixels with indeterminate a c and 
T < T t — 3 K are hereafter referred to as “dark” pixels. 
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282 281 280 279 

1 1 .5 /xm TEMPERATURE (K) 

Fig. 4a. VIS-IR histogram of GOES pixels over FMC at 1500 
UTC 27 October 1986 ( numbers denote frequency of occurrence of 
temperature-count pairs). 


They are not used to sol ve ( 7 ) . Their impact and origins 
are discussed in section 4d. The cloud emittances 
are plotted in Fig. 5a against the measured reflectances 
for the case in Fig. 4b. Eliminating the dark pixels and 
applying (5) to the data in Fig. 5a yields the cloud 
albedo values. Averaging the emittances for a given 
albedo produces the mean values and the standard de- 
viations plotted in Fig. 5b. The solid line represents 
the solution to (7) using the average value of k = 5.1. 
The mean beam emittance and VIS optical depths for 
the data used in Fig. 5b are 7 b = 0.38 and r v = 0.59. 

The data from FMC in Fig. 5b for this hour are 
compared to those from MAD and WAU shown in 


Fig. 6. Apparently, the clouds over WAU are much 
denser than those over MAD, while the MAD obser- 
vations are similar to those over FMC. Values of T c 
differed by only 1 K among the sites, while T t ranged 
from 225 K at WAU to 2 1 7 K at MAD. Depolarization 
ratios derived from the lidar returns indicated inter- 
mittent liquid layers during the morning of 28 October, 
especially at ~ 1 500 UTC. Those liquid layers may be 
the source of the larger emittances over WAU. Data 
from all three sites were combined, averaged, and fit 
with (7) yielding k = 5.6. The scatter in the means 
between the sites at a given hour is of the same order 
as that for different hours at the same site as seen in 
Fig. 7 for FMC at 3 times during 28 October. 

1 ) GOES-surface lidar results 

A summary of the results for the case study 27-28 
October is given in Table 2. Cloud-top heights range 
from 9.5 to 1 1 .0 km at all three sites. The cloud-center 
temperatures vary by about 25 K. Cloud optical depths 
were much greater over WAU than over the other sites. 
Dark pixels were found more often over FMC and 
MAD' than over WAU. 

Due to dropouts, the only data available for 1 500 
UTC during the primary IFO cirrus days occurred on 
28 October. Data from other days were available for 
most of the afternoon hours. The combined datasets 
permitted coverage of the full range of emittances at a 
given hour as illustrated in Fig. 8 for 2000 UTC. Data 
for all of the hours used in the case study and the entire 
IFO analyses are shown in Figs. 9a and 9b, respectively. 
The clouds over the area during 27-28 October were 
generally thinner with lower emittances than most of 
those observed during the other days. The combined 
datasets (Fig. 9b) yield a large number of samples for 
e b < 0.5 and « b > 0.8 and relatively few for intermediate 
values of f b . Apparently, the cirrus, which occurred 
during the IFO, tended to be either very thick or rel- 
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Fig. 4b. VIS-IR histogram of GOES pixels over FMC at 1500 UTC 28 October 1986. 
Dashed line refers to maximum clear-sky count. 
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Fta. 5a. Cloud emittances and observed reflectances for T < T, Fig. 6. Average doud emittance versus cloud albedo at 1500 UTC 
- 3 K derived from Fig. 4b. Dashed line refers to maximum clear- 28 October 1986 over three sites from GOES VIS-IR data, 
sky reflectance. 


atively thin. Perhaps, a longer time period would pro- 
duce more uniform sampling. No data are found for 
tb < 0.08 since no clouds are recognized if T > T s 
-3 K. 

Table 3 summarizes the values of r v derived from 
all of the GOES-surface IFO data (Fig. 9b) and from 
case study data only (Fig. 9a) for each relevant time. 
The average scattering angles, 0, between the sun, sat- 



FiG- 5b. Cloud albedos and mean cloud emittances derived from 
Fig. 5a without “dark” pixels. Vertical lines represent standard de- 
viations. The curved line denotes the regression fit to Eq. (7). 


ellite, and scene are also listed in Table 3. Visible optical 
depths observed during the case study are less than half 
of those observed for all of the IFO cirrus days. 

The mean vertical emittance is given as a function 
of temperature in Fig. 10. The dashed line correspond- 
ing to the results from Fig. 7a of Platt et al. ( 1987) is 
included for comparison. Although there is a general 
increase in t a with increasing cloud temperature, the 



Fig. 7. Average cloud emittance and albedo at three times during 
28 October 1986 over FMC from GOES VIS-IR data. 
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Table 2. Observed and computed cloud properties for October 1986 case study. 



Time 

T, 

z x 

h 

T' 

T, 



Dark pixel 


t 

Site (Day) 

(UTC) 

(K) 

(km) 

(km) 

(K) 

(K) 



(%) 

Ty 

FMC (27) 

2030 

287.7 

9.6 

1.1 

227.8 

225.1 

0.11 

0.23 

4 

0.20 

U 

2100 

286.9 

10.3 

1.6 

228.7 

220.3 

0.19 

0.52 

0 

0.55 

2X 

(28) 

1330 

278.3 

9.5 

2.1 

237.9 

227.7 

0.22 

0.12 

0 

0.22 

l.( 

1400 

279.7 

11.0 

1.3 

231.6 

216.8 

0.14 

0.16 

5 

0.21 

U 


1430 

280.5 

11.2 

1.5 

230.0 

215.6 

0.18 

0.32 

11 

0.38 

2.( 


1500 

281.6 

10.1 

1.3 

239.4 

222.5 

0.31 

0.62 

7 

0.65 

l.< 


1600 

283.0 

10.0 

3.0 

240.1 

223.2 

0.80 

1.77 

0 

1.61 

2.1 


1700 

287.9 

10.4 

1.9 

223.0 

219.6 ' 

0.23 

0.95 

0 

0.74 

3 X 


1900 

285.4 

10.8 

2.1 

227.5 

217.1 

0.23 

0.64 

0 

0.45 

1.? 


1930 

285.5 

10.9 

2.6 

228.9 

216.6 

0.41 

0.99 

3 

0.79 

l.< 


2000 

285.9 

11.0 

1.9 

224.4 

216.3 

0.41 

1.05 

0 

0.87 

l.< 


2030 

284.8 

10.9 

1.7 

220.7 

216.4 

0.62 

1.78 

6 

1.55 

l.< 


2100 

280.9 

10.7 

2.4 

230.7 

216.7 

0.79 

1.41 

5 

1.53 

2.< 


2130 

281.2 

10.8 

2.2 

226.3 

216.4 

0.44 

0.75 

11 

1.33 

i: 


2200 

276.8 

10.9 

2.4 

221.7 

216.1 

0.32 

0.30 

0 

0.87 

2j 

MAD (28) 

1330 

278.4 

■a 

3.4 

237.9 

217.9 



3 


2.: 

1500 

280.2 

10.8 

4.2 

242.8 

217.7 

0.24 

0.65 

2 


2X 


1600 

281.7 

10.6 

3.7 

241.9 

218.5 

0.58 

1.18 

0 

1.07 

l.< 


1700 

284.8 

10.5 

4.3 

240.4 

218.8 

0.40 

0.98 

0 

0.75 

l.< 


1800 

289.1 

9.8 

3.2 

229.4 

224.6 

0.14 

0.95 

0 

0.75 

5.: 


1930 

287.7 

10.7 

1.8 

222.1 

217.3 

0.10 

0.26 

3 

0.13 

i.: 


2000 

286.5 

10.7 

1.8 

221.7 

217.1 

0.10 

0.55 

7 

0.43 

4/ 


2030 

286.6 

10.7 

3.1 

223.3 

217.1 

0.34 

0.86 

3 

0.76 

2.; 


2100 

286.0 

10.5 

2.4 

223.3 

217.2 

0.29 

0.56 

33 

0.59 

i.< 


2130 

284.5 

10.8 

4.0 

232.2 

216.4 

0.55 

0.46 

3 

0.73 

i.: 


2200 

281.1 

10.5 

4.2 

237.4 

216.8 

0.44 

0.27 

13 

1.24 

2. 

WAU (28) 

1500 

279.0 

9.8 

3.2 

238.0 

225.0 

1.67 

3.20 

0 

3.50 

2.( 

1600 

282.5 

10.1 

3.1 

238.0 

222.4 

1.40 

2.42 

0 

2.19 

1.: 


1700 

287.1 

10.1 

3.9 

235.0 

222.3 

1.21 

2.52 

0 

1.91 

1.' 


1800 

288.4 

10.3 

4.1 

231.0 

220.3 

1.13 

3.13 

0 

2.35 

4.( 


1900 

285.7 

11.0 

0.8 

217.1 

216.4 

0.19 

0.77 

0 

0.55 

3. 


1930 

283.9 

10.9 

2.1 

226.4 

216.6 

0.30 

1.03 

0 

0.82 

2.* 


2000 

284.4 

10.9 

2.9 

236.1 

216.5 

1.04 

2.42 

0 

1.93 

i.; 


2030 

285.4 

10.7 

3.1 

231.0 

216.9 

1.34 

3.73 

0 

3.11 

2.1 


2100 

284.1 

10.6 

3.2 

234.0 

217.0 

0.36 

1.43 

7 

1.66 

3.; 


2130 

279.7 

11.0 

1.0 

217.0 

216.0 

0.29 

0.45 

1 

0.76 

2. 


values derived for the case study are generally lower 
than the IFO means. The case study averages parallel 
those found by Platt et al. ( 1987) using ground-based 
observations in Australia. The mean value of e, found 
here for the case study results is 0.33, the same as that 
from Platt et al. ( 1987). 

Goud thickness versus T c is shown in Fig. 1 1 . Max- 
imum average thickness is found at T c = 232 K for 
the IFO and at T c »» 241 K for the case study. The 
case study values of h increase nearly monotonically 
with T c . The IFO results show a tendency towardis 
thinner clouds for T c < 220 K and for T c > 240 K. 
The thickest clouds were observed during 22 October 
and 1 November. During the case study, the observed 
clouds were 2.6 km thick on average compared to a 
mean thickness of 4.2 km during the remaining cirrus 
days of the IFO. 

Goud IR volume absorption coefficient, o, = r t /h, 
is plotted against T c in Fig. 12. The results reveal a 
trend of increasing extinction with T c in the IFO data 


(correlation coefficient of 0.71 ). The average value < 
o e is 0.20 for the case study data. 

The variation of vertical emittance with T c for tl 
IFO data differs from that found for the case study ai 
Platt et al. ( 1987), probably because of substanti 
differences in cloud depths. This difference in thickne 
is evident in Fig. 1 1 where the case study data, a subs 
of the IFO data, are also consistent with the results 
Platt et al. (1987), probably because of substanti 
ume absorption coefficients in Fig. 1 2 are not as w< 
behaved as emittance and cloud thickness when co 
sidered as a function of T c . Except for the highest ar 
lowest temperatures (which are poorly sampled), 
for the case . udy data is close to that observed by Pis 
et al. (1987). Differences between the IFO and ca 
study data are also seen in the behavior or o t suggest i> 
that thickness is not the only discrepancy between tl 
case study and other IFO clouds. Vertical distribution 
of nongaseous cloud water content and particle shapt 
sizes, and phase may also affect these differences. 
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Fig. 8. Average cloud emittances and albedos from GOES at 2100 
UTC over FMC during different IFO days. 


FIG. 9b. Same as Fig. 9a, except for all IFO scenes. 


2) GOES/ER-2 cloud parameters 

Table 4 lists the data and derived parameter values 
for the five cirrus matches between the GOES and ER- 
2 lidar for the case study days. Data taken during 27 
October were located just outside the northwestern 
boundary of the IFO area. On 28 October the ER-2 
was over Lake Michigan at ~ 1600 UTC, while it was 
near FMC at 1 800 UTC. The three scenes at 1 600 UTC 
correspond to small contiguous areas. The small optical 



Flo. 9a. Average cloud emittance versus cloud albedo for all case 
study scenes from GOES over surface sites where T, = T c . 


depths indicate that the aircraft also sampled relatively 
thin clouds during the case study. 

3) Cloud parameters from AVHRR data 

Figure 13 presents the averaged cloud emittance- 
albedo pairs derived from both AVHRR and GOES 
data taken for the 0.25° X 0.35° area over WAU on 
28 October. Although the relative viewing angles were 
within ±8° for the two satellites, the absolute viewing 
conditions were different. The GOES viewed the region 
from the south, while the AVHRR viewed from the 
north with the solar plane almost midway between the 
two views. Discrepancies in the albedo range may be 
due to this misalignment (some dark pixels were found 
in the GOES results). Resolution differences, however, 
would tend to produce a smaller range in both albedo 
and emittance for the GOES (8 km) relative to the 
AVHRR ( 1 X 4 km 2 ) radiances as observed in Fig. 
13. In general, the GOES data are very close to the 
AVHRR results except for the smaller values of a c . 
The latter may be affected by the presence of “dark” 
pixels. The values for the AVHRR VIS and IR optical 
depths are ~ 10% greater than their GOES counter- 
parts. 

Another comparison of GOES and AVHRR results 
is given in Fig. 14 for data taken over WAU 2 Novem- 
ber. The AVHRR viewed the scene with 8 = 18° and 
i ' = 38°. The squares correspond to AVHRR data as 
reduced for the given viewing angles. AVHRR emit- 
tances corrected to the GOES viewing zenith angle us- 
ing (2) are denoted with the crosses. The range in a c 
is smaller for GOES than for the AVHRR, consistent 
with the resolution differences. In this case, the 
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Table 3. Reflectance parameters computed for all GOES-surface lidar data. 








Nominal (using X c ) 



Reanalyzed (using X’ € ) 






e 

All data (IFO) 

Case study 

All data (IFO) 

Case study 

i i me 

Cases 

8 










(UTC) 

(IFO) 

n 

n 

n 

Ty 

t 

Ty 

« 

Ty 

{' 

Ty 

«• 

1330 

4 

80.2 

106.6 

109 

1.28 

1.17 

0.15 

0.83 

0.28 

1.61 

0.28 

1.6 

1400 

3 

75.6 

112.6 

117 

2.21 

1.42 

0.17 

1.12 

3.22 

2.32 

0.21 

1.5 

1430 

2 

71.3 

118.6 

124 

1.07 

1.79 

0.32 

1.72 

1.35 

2.22 

0.38 

2.0 

1500 

3 

67.9 

125.7 

131 

1.72 

2.02 

1.72 

2.02 

1.87 

2.15 

1.87 

2.1 

1600 

5 

61.0 

140.0 

146 

2.93 

2.50 

1.82 

2.11 

2.64 

2.25 

1.64 

1.9 

1630 

2 

57.8 

147.7 

153 

2.59 

1.85 





3,06 

2.18 



1700 

6 

57.6 

156.6 

160 

2.58 

2.90 

1.57 

2.93 

1.96 

2.21 

1.20 

2.2 

1800 

6 

57.3 

174.1 

173 

2.43 

2.96 

1.95 

6.21 

1.92 

2.33 

1.48 

4.8* 

1830 

1 

56.1 

177.1 

175 

3.65 

2.06 

— 



3.88 

2.19 



1900 

7 

59.9 

169.0 

168 

1.62 

3.12 

0.69 

3.43 

1.82 

2.12 

0.49 

2.3< 

1930 

4 

61.1 

160.8 

162 

1.73 

2.74 

0.75 

2.76 

1.34 

1.97 

0.57 

1.9 

2000 

8 

65.0 

153.6 

154 

2.07 

2.68 

1.69 

2.87 

2.27 

2.26 

1.35 

2.31 

2030 

5 

68.3 

146.5 

147 

1.56 

2.55 

1.87 

2.37 

1.33 

2.23 

1.59 

2.0: 

2100 

6 

72.8 

140.3 

140 

1.07 

2.07 

1.02 

2.19 

1.17 

2.24 

1.12 

2.31 

2130 

4 

76.2 

133.6 

. 132 

1.12 

1.29 

0.55 

1.25 

0.81 

2.23 

0.94 

2.0) 

2200 

5 

81.2 

128.1 

125 

0.80 

0.84 

0.28 

0.70 

1.04 

2.30 

0.80 

2.3< 

Totals and 













means 

71 

67.2 

145.4 

— 

1.69 

2.13 

1.04 

2.17 

1.80 

2.19 

1.11 

2.K 


AVHRR data produce a much lower minimum cloud 
albedo. The application of (2) appears to have pro- 
duced very similar emittances for the two datasets, al- 
though there is a 20% difference in the average values 
of r e . The mean VIS optical depths differ by a factor 
of 2. 

All of the coincident AVHRR and GOES data are 
summarized in Table 5. The AVHRR IR optical depths 



FIG. 10. Variation of cloud vertical emittance with cloud-center 
temperature for all IFO (circles) and case study (squares) data from 
GOES over surface lidar sites. Vertical lines denote standard devia- 
tions. Dashed line adapted from Platt et al. (1987). 


are consistently greater than the corresponding GOES 
values by ~0. 1. This difference indicates the possibility 
of a calibration offset in the thermal channels. Despite 
this obvious bias, the good relative agreement in Fig. 
14 between the corrected AVHRR emittances and the 
GOES emittances suggests that ( 3 ) is a reasonable ap- 
proximation to the IR absorption optical depth. Any 
IR scattering effects that are ignored here are apparently 



FIg .11. Variation of cloud thickness with cloud-center temperature 
for all IFO (circles) and case study (squares) data from GOES over 
surface lidar sites. Vertical lines denote standard deviations. 
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CLOUD TEMPERATURE (K) 


Rg . 1 2. Variation of infrared volume extinction ( absorption ) coef- 
ficient with cloud-center temperature for all IFO (circles) and case 
study (squares) data from GOES over surface lidar sites. Vertical 
lines denote standard deviations. Dashed curve adapted from Platt 
etal. (1987). 


insignificant compared to the other error sources. Dif- 
ferences between the GOES and AVHRR values of r v 
vary from scene to scene. Even when the times and 
angles are very close and the data appear similar as in 
Fig. 13, there are substantial differences in r v . Over 
FMC during 28 October, there is good agreement be- 
tween the parameters however. The outstanding dif- 
ferences may be attributable to a number of factors 
that are discussed in section 3e. 

d. Results and discussion for cloud-top temperatures 
with GOES 

The analyses discussed above were also performed 
for the GOES-derived emittances for T z = T t . A plot 
of all of the mean cloud emittance-albedo pairs is 
shown in Fig. 15. The largest concentrations of data 
are found for t b ( T t ) < 0.8. In general, a c is greater for 
a given value of < b than it is in Fig. 9b. Goud-top emit- 


tance tends to plateau at < b = 0.86, while cloud-center 
emittance appears to level at t b 0.98. The lower 
emittances lead to diminished values of r e relative to 
those derived for T c . 

The emittance ratio, r, = « b ( T t )/( b ( T c ), was com- 
puted for discrete intervals of T c . Mean values and 
standard deviations of these ratios are shown in Fig. 
1 6. The emittance ratio increases almost linearly with 
decreasing cloud center temperature. Standard devia- 
tions about a given mean ratio are less than 0. 1 . The 
emittance ratio is close to unity for T c < 2 1 5 K. 

The well-correlated variation of r, suggests the pos- 
sibility that T x as well as T c may be retrieved with VIS- 
IR radiance pairs. This ratio integrates many of the 
other parameters examined earlier. For example, cloud 
thickness in Fig. 1 1 is least at the highest altitudes and 
increases before leveling or even decreasing for tem- 
peratures around 235 K. For the highest clouds, there 
is little difference between T c and T t because the clouds 
are not very thick. Since cloud depths are greater at 
lower altitudes, it is possible to sense radiation from 
areas deep within the cloud thereby causing greater 
differences between T c and T t . As the depth of the cloud 
decreases, the ratio should approach unity. At higher 
temperatures (T c > 250 K), the relationship of r, to 
T c may not be as well defined because liquid water 
becomes more common and the mean cloud depth 
may not be dependent on T c . Whether the relationship 
shown in Fig. 16 is typical for all observing angles is 
also unknown. Additional sampling from other angles 
and over a wider variety of temperatures would help 
to better define the relationship between r, and T c . 
However, the results shown in Figs. 11, 12, and 16 
suggest that it may be possible to obtain reasonable 
estimates of h and T, over a limited range of T c . 

e. Emittance uncertainties 

The parameter values derived here are subject to 
considerable uncertainty as evidenced by the results in 
Table 5 and the large standard deviations in earlier 
figures. Potential sources of error abound in an analysis 
of this type due to the large number of variables and 
the nonuniformity of cirrus clouds. 

Parameters derived from the lidar essentially provide 


Table 4. Observed and computed cloud properties for case study ER-2 data. 


Day 

Time 

(UTQ 

Lat. 

(°N) 

Lon. 

(*W) 

T. 

(K) 

T c 

(K) 

h 

(km) 

u 

T v 

t 

Dark 

pixel 

(%) 

Tv 

r 

27 * 

1830 

45.8 

93.1 

290.0 

225.0 

0.5 

0.12 

0.69 

5.73 

0 

0.46 

3.86 


1900 

45.3 

92.5 

288.5 

234.0 

0.5 

0.14 

0.51 

3.43 

0 

0.39 

2.64 


1930 

44.9 

91.1 

288.4 

234.0 

0.5 

0.24 

0.73 

2.91 

3 

0.61 

1.68 

28 

1600 

44.6 

87.0 

283.0 

229.0 

3.7 

0.43 

0.86 

1.99 

0 

0.84 

1.95 


1600 

44.5 

87.1 

283.0 

229.0 

3.7 

0.42 

0.92 

2.12 

0 

0.91 

2.10 


1600 

44.5 

87.1 

283.0 

229.0 

3.7 

0.45 

1.00 

2.21 

0 

0.99 

2.21 


1800 

43.6 

89.4 

288.9 

228.0 

1.5 

0.28 

0.90 

3.23 

0 

0.69 

2.57 
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Fig. 13. Comparison of cloud albedos and emittances derived from 
28 October 1986 GOES and AVHRR data taken at ~2100 UTC 
over WAU. 


a two-dimensional view of the cirrus clouds. The va- 
lidity of the assumption that z c ,h, and z, represent the 
average cloud heights within the large areas covered by 
the strip of pixels is difficult to evaluate. One means 
of estimating how well the lidar data represent the large- 
area cloud characteristics is to examine the differences 
between the strip of pixels and surrounding areas. The 
rms difference between the emittances for the strip and 
the box containing the strip is 0.05 or 7%. This differ- 
ence is equivalent to a ±0.7 km variation in cloud- 
center height between the strip and the box. Changes 
of 2 km in cloud-center altitude during a given half 
hour are common as seen in Fig. 1 . The variations in 
the small-scale lidar data are greater than those in the 
large-scale satellite data as expected. 

Assuming that the large-scale differences are repre- 
sentative of the lidar-satellite scale differences, it is es- 
timated that the use of lidar data to set z c causes an 
uncertainty in of ±10% based on an average value 
for t, of 0.62. Note that no clouds with « b < 0. 1 were 
included in the analysis because of the cloud-detection 
threshold of 3 K. A conservative estimate of the un- 
certainty in T t is ±2 K. Inclusion of this error raises 
the overall uncertainty in e, to ± 1 3%. This uncertainty 
in is equivalent to a ±20% uncertainty in r e for a 
given scene over the range of optical depths considered 
here. The AVHRR-GOES comparisons are, on aver- 
age, within this uncertainty level. The average IR op- 
tical thickness is 0.96 for all 7 1 scenes. From the strip 
and box comparison, it is also estimated that z, and h 
have uncertainties of ±0.7 km. 

Another source of uncertainty in t, is the use of a 
mean cloud height for the entire scene. This error 
source may be examined by performing a pixel-by-pixel 
analysis on a scene that varies systematically with time. 
One example is the cloud over FMC between 2020 and 


2050 UTC. The GOES pixels from the corresponding 
wind strip data were averaged in lines perpendicular 
to the wind vector. Using the wind speed, these aver- 
aged pixels were converted to times and aligned with 
the lidar-defined cloud parameters. The results shown 
in Fig. 17 indicate good alignment between the two 
datasets. In this case, it appears that the lidar data pro- 
vide an accurate cross section of the cloud. The GOES 
reflectance increases as the cloud thickens and T in- 
creases as z c lowers. Equations ( 1 ) and ( 5 ) were applied 
to each average pixel using T c derived from Fig. 17 to 
determine t c and r v . Figure 18 shows the variation of 
the parameters with time. Although the thin part of 
the cloud is detected with the IR data, a value for r v is 
not computed since the reflectances are lower than that 
for clear skies. Nevertheless, the mean values for r e 
and <b derived on a pixel-by-pixel basis are 0.59 and 
0.54, respectively, compared to 0.62 and 0.59 derived 
for the entire scene using a mean value of z c . This 
comparison suggests that the error in e b for using the 
mean cloud height is around 10% with a slight tendency 
to bias the values to the high end because of nonlinear 
effects. While these results may not represent all cases, 
they indicate that the use of a mean cloud height for 
the analysis is a reasonable approach. 

4. Relationship between VIS and IR parameters 

a. Scattering efficiency ratio 

For a given cloud particle with cross-sectional area 
lira 2 , the VIS scattering cross section is 

j8v = QsoZira 2 , 

and the IR absorption cross section is 
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Fig. 14. Same as Fig. 13, except for 
2 November 1986 at ~2000 UTC. 
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Table 5. Comparison of AVHRR- and GOES-derived cloud parameters. 


AVHRR Parameters 


GOES Parameters 


Site 

(Day) 

T. 

(K) 

Time 

(UTC) 

e 

n 

u 

Ty 

{ 

r 

Dark 

pixel 

% 

Time 

(UTC) 

r e 

r v 

f 

Dark 

pixel 

% 

Ty 

£' 

FMC 

(28) 

231 

2100 

141 

0.83 

1.47 

1.95 

1.96 

0 

2100 

0.72 

1.36 

1.79 

3 

1.48 

1.94 

(30) 

223 

2042 

127 

0.25 

0.35 

1.13 

1.79 

18 

2030 

0.17 

0.63 

3.50 

3 

0.57 

3.12 

(2) 

230 

2006 

92 

0.49 

1.18 

1.94 

1.36 

5 

2000 

0.39 

0.97 

:.42 

0 

0.80 

1.98 

WAU 

(28) 

234 

2100 

142 

0.31 

1.69 

3.50 

21.4 

0 

2100 

0.26 

0.78 

2.20 

6 

0.89 

2.41 

(2) 

230 

2006 

98 

1.34 

4.61 

3.21 

1.78 

0 

2000 

1.16 

3.26 

2.72 

0 

2.48 

2.09 


& = C?abs2ira 2 , 

where the VIS scattering efficiency is Q the IR ab- 
sorption efficiency is <2«ts> and the effective particle 
radius is a. The extinction optical depth for a given 
wavelength, A, is 

r(X) = /S(A) j N(z)dz, 

where the integral is over some finite height, 0 is the 
extinction optical depth, and N(z) is the number den- 
sity of particles at z. If it is assumed for a given cloud 
that the VIS extinction is due entirely to scattering, the 
IR extinction is due entirely to absorption, and the 
extinction coefficients are constant, then the scattering- 
to-absorption ratio is 

£ = (?sca/(2abs = |8v/& = T v /T e . 

Therefore, if £ and r v are known, then the value of r e 
may be determined. This efficiency ratio forms the basis 



of the ISCCP cloud-height correction algorithm. The 
values of O^, and Q»bs depend primarily on the particle 
size in terms of the size parameter, 2 ira/A, and the 
index of refraction. Mie scattering calculations of Q ^ 
and Cabs indicate that £ for the wavelengths of interest 
may vary from ~5 for a = 2.0 //m to ~2.4 for a = 10.0 
(im (Sassen 1981 ). The value of £ asymptotically ap- 
proaches 2.0 for very large particles; i.e., 2ir a/X > 1 
(e.g., Hansen and Travis 1974). Platt et al. (1987) have 
presented some theoretical results that indicate £ is also 
a function of particle shape and may have values 
slightly less than 2.0. 

For each satellite-lidar dataset in this study, the scat- 
tering efficiency ratio is computed for all of the pixels, 
N, having a discrete value of p as follows. 

N 

£ - 2 (tv/tJ/N. ( 8 ) 

i-1 



CLOUD TEMPERATURE (K) 

Fig. 16. Variation of ratio of cloud-top emittance to cloud-center 
emittance on cloud-center temperature for all IFO (circles) and case 
study (squares) data from GOES over surface lidar sites. Vertical 
lines denote standard deviations. 
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Flo. 17. Comparison of GOES and lidar observations along wind 
vector over FMC during 28 October 1986. 


Only one temperature, T c , is used to compute r e for a 
given dataset since only one average cloud height is 
derived for each time. Changes in the actual cloud 
height and thickness within the scene (e.g., Fig. 1 ) tend 
to introduce variations in r e for a given reflectance. 
Thus, the mean value of £ is computed for each cloud 
reflectance value to minimize the effects of cloud height 
variability. 

Visible optical depths for the case study GOES data 
are given in Table 2. Tables 3-5 summarize the results 
of applying ( 8 ) to all of the data. Despite the differences 
between the values of r v for IFO and case study results 
in Table 3, two similarities are quite evident in a com- 
parison of the respective scattering efficiencies. For both 
datasets, the scattering efficiency appears to increase 
with decreasing 0 o and increasing 8. At high values of 
Oo, £ is well below the expected limit of 2. The average 
values of £ are also very close, 2.17 and 2.13, for the 
case study and IFO, respectively. In Table 4, the greatest 
values of £ also occur near local noon ( 1800 UTC). 

b. Reanalyzed visible data 

The temporal dependencies of r v and £ are not re- 
alistic. They are primarily due to shortcomings in the 
analysis treatment of Xc. The anisotropy of the reflected 
radiation field for real clouds depends on the optical 
thickness, incident radiation, microphysical properties 
of the cloud, and the morphology of the cloud field. 
The value of X c used here is fixed for a given set of 
angles and represents an empirical average for all cloud 
types. The average cloud optical depth in the bidirec- 


tional reflectance model used here is probably close 
10, while t v for the clouds analyzed in this study 
generally smaller than 2. Since cirrus clouds are t 
only type considered here and the angles are fixed f 
a given hour, it is likely that Xc will be biased wi 
respect to local time. There will also be random errc 
in x c due to variations in microphysics, morpholog 
and cloud optical depth for a given hour. The ma 
nitudes of these errors are currently unknown, but a 
potentially large. Assuming that the time samplii 
represents a random sampling of X c , the averages 
various parameters derived from all times should 1 
relatively unbiased. 

Using the assumption that the mean value of £ 
independent of time, new values of Xc were determim 
from (5) using the observations of r e and the mea 
scattering ratio of 2. 1 3. These new values were average 
at each time. The means, denoted as Xc, are given i 
Table 6 with the mean nominal values from the bid 
rectional reflectance model. The results indicate moi 
anisotropy in the cirrus reflectance pattern than in tl 
empirical model. The data of Takano and Liou ( 19891 
indicate that reflectance anisotropy diminishes wit 
increasing optical depth for scattering by hexagon; 
columns. Thus, the smaller optical depths of the clouc 
here compared to those for the clouds in the empiric; 
model are probably responsible for the larger range c 
Xc compared to Xc- 

To eliminate the temporal (angular) dependence c 
r v , the data were reanalyzed using Xc. The values c 
the visible optical depth and the scattering ratios dt 
termined from the reanalysis are listed in Tables 2- 
and denoted with the primed variables, t' v and £', re 
spectively. Values of £' for the AVHRR data wer 
computed by changing the nominal values of X c so tha 
the derived value of t' v for the AVHRR data was equa 
to t'v for the corresponding GOES data. In general, th 
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Fig. 18 . Cloud optical properties derived from Fig. 17 . 
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Table 6. Anisotropic reflectance factor comparison. 


Time 

Number 

Xc 

x; 

Mean 

Standard 

(UTC) 

of pixels 

(nominal) 

({ = 2.13) 

difference 

deviation 

1330 

213 

0.823 

0.618 

0.205 

0.104 

1400 

177 

0.864 

0.735 

0.129 

0.050 

1430 

113 

0.896 

0.809 

0.087 

0.113 

1500 

177 

0.922 

0.896 

0.026 

0.079 

1600 

215 

0.980 

1.034 

-0.054 

0.123 

1630 

146 

1.002 

0.912 

0.090 

0.044 

1700 

212 

1.020 

1.195 

-0.175 

0,264 

1800 

158 

1.050 

1.206 

-0.156 

0.348 

1830 

72 

1.061 

1.025 

0.036 

0.000 

1900 

262 

1.029 

1.214 

-0.185 

0.220 

1930 

236 

1.008 

1.147 

-0.139 

0.206 

2000 

245 

0.982 

1.091 

-0.109 

0.164 

2030 

290 

0.953 

1.008 

-0.055 

0.103 

2100 

222 

0.918 

0.890 

0.028 

0.093 

2130 

283 

0.882 

0.713 

0.069 

0.154 

2200 

210 

0.843 

0.579 

0.264 

0.082 


new mean optical depths in Tables 2 and 3 have been 
increased for hours near the terminator and decreased 
near noon. The temporal dependence of the scattering 
efficiency ratios has been effectively eliminated for the 
IFO averages, while some tendency for maximum val- 
ues near noon remains for the case study results. Re- 
sulting means of £' are not equal to 2.13 because of 
the nonlinearities in the solution of ( 5 ) and ( 8 ) and 
because four additional scenes were eliminated due to 
filtering. In Table 4, all of the optical depths and ratios 
were decreased yielding a mean value for {' of 2.43. 
The new values of X* for the AVHRR data (Table 5) 
reduced the differences between the GOES and 
AVHRR scattering efficiency ratios with respective 
means of 2.3 1 and 2.35 for the five values of £'. The 
greatest difference between the two datasets is found 
for FMC during 30 October when almost 20% of the 
AVHRR data are classified as dark pixels. 

Mean values of r' v were computed for discrete in- 
tervals of cloud thickness. The results shown in Fig. 
19 reveal a strong, almost linear, dependence of t' v on 
cloud thickness. This apparent relationship lends fur- 
ther support to the potential for deriving cirrus depth 
from the VIS-IR radiance data. Figure 20 shows the 
variation of mean £' with IR optical depth. The stan- 
dard deviations of the IFO data are denoted with the 
vertical lines. The case study averages (squares) vary 
from 1.79 to 2.54, while the IFO means (circles) range 
from 1.91 to 2.33. There is good consistency between 
the IFO and case study results for r e < 1 .0. Gouds with 
r e > 1.0 were undersampled during the case study (Fig. 
9a). These results show no indication of a dependence 
of ( ' on t c . 

Average scattering efficiencies were also computed 
for discrete levels of T c (Fig. 21 ). For T c > 230 K, £' 

2.0. Despite the large standard deviations, the mean 
value of £' appears to increase with decreasing tem- 
perature for T c < 230 K. At T c = 217 K, the mean 



Fig. 19. Variation of cloud visible optical depth with cloud thick- 
ness for all IFO (circles) and case study (squares) data from GOES 
over surface lidar sites. Vertical lines denote standard deviations. 


scattering efficiency ratio reaches a maximum of ~2.6 
for the IFO data. To determine if there is a significant 
difference between the results for cold and warm cirrus 
clouds, average values for £' were computed for T c 
230 K and for T c > 230 K. For the thirty cases with 
T c < 230 K (cold clouds), the mean value of £' is 
2.43 with a standard deviation of 0.89. For the other 



Fig. 20. Variation of mean scattering efficiency ratio with infrared 
optical depth for all IFO (circles) and case study (squares) data from 
GOES over surface lidar sites. Standard deviations denoted by vertical 
lines. 
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37 scenes (warm clouds), the mean scattering efficiency 
ratio is 2. 1 1 with a standard deviation of 0.7 1 . Appli- 
cation of sampling theory with an assumption of nor- 
mal distributions yields £' = 2.43 ±0.17 for the cold 
clouds and £' = 2.11 ± 0. 1 2 for the warm clouds. ( The 
mean for all 67 cases, 2.25, differs from that in Table 
3 because of different averaging techniques.) Using a 
two-tailed Student /-test, it was found that the cold- 
cloud (T c ^ 230 K) ratios are different from the warm- 
cloud values at the 90% confidence level. The mean 
ratio for the warm clouds is not significantly different 
from 2.0, the lower limit expected for large particles. 

As noted earlier, the value of £ is highly dependent 
on the particle size. Neglecting the effects of particle 
shape, it may be concluded from these results that large 
particles are the dominant scatterers in the clouds for 
T c > 230 K. At colder cloud temperatures, smaller 
particles (e.g., a < 20 pm ) apparently contribute more 
to the scattering. Thus, the average particle size in the 
colder clouds must be significantly less than the warmer 
clouds. Figure 2 1 suggests that the proportion of small 
particles in the clouds increases as the cloud-center 
temperature decreases. These findings are consistent 
with the results of Platt and Dilley (1981) who deter- 
mined that the lidar-measured backscatter-to-extinc- 
tion ratio is relatively constant for T c > 233 K. The 
value of that ratio changes abruptly to a new value at 
temperatures slightly less than 233 K. Heymsfield and 
Platt ( 1984) were able to show that the number density 
of large particles in cirrus clouds decreased substantially 
when the temperature dropped below -40°C. The VIS 
extinction coefficients derived from the GOES data are 



Fig. 21. Variation of mean scattering efficiency ratio with cloud- 
center temperature for all IFO ( circles ) and case study ( squares) data 
from GOES over surface lidar sites. Vertical lines denote standard 
deviations. 


0.41 km "‘ and 0.45 km for T c « 230 K and T c > 230 
K, respectively. These values fall in the range of values 
computed by Heymsfield and Platt ( 1984), but show 
less dependence on temperature. 

Although particle size is very important in the de- 
termination of £, Platt et al. (1987) have presented 
theoretical results indicating that £ is also a function 
of particle shape. The results of Heymsfield and Platt 
( 1984) show that ice crystal habits change significantly 
when the cloud temperature drops below 233 K. Thus, 
both factors, smaller crystals and shape changes, may 
explain the increase in £ for T c < 230 K observed in 
the current dataset. Microphysical measurements using 
probes sensitive to small ice particles are needed to 
further clarify these results. 

c. Assumptions for visible optical depth calculations 

The determination of Tv relies on the assumptions 
that the clouds can be modeled as plane parallel entities, 
that (5) is an adequate parameterization of a more 
exact radiative transfer model, and that the cloud al- 
bedo model, or the implied scattering phase function, 
is appropriate for the clouds studied here. The lidar 
plots in Figs. 1 and 2 show considerable structure in 
the clouds observed during 28 October. Such mor- 
phology is probably typical of most of the cirrus clouds 
observed during the IFO. Thus, the clouds do not ap- 
pear to be plane-parallel and homogeneous. Because 
exact modeling of the nonhomogeneous clouds is not 
possible at this time, no attempt is made to estimate 
the effects of cloud morphology on the derived values 
of VIS optical depth. 

One means of evaluating the parameterization is to 
compare the results of ( 5 ) with values of t v derived for 
each scene with a detailed radiative transfer model us- 
ing the scattering phase function that produced the 
cloud albedo model applied in ( 5 ) . Takano et al. ( 1 989 ) 
computed t v for several of the IFO scenes using a three- 
level adding-doubling radiative transfer model with 
two different scattering phase functions. One of the 
phase functions was the same as that used to derive 
the cloud albedo model employed here. The other was 
for ice spheres with areas equivalent to the hexagonal 
columns used in the first phase function. Takano et al. 

( 1989) used a Lambertian surface with an albedo of 
0. 1 1 to simulate the Wisconsin land areas. They de- 
termined the value of r» for 16 scenes during the case 
study period by matching the reflectance of the adding- 
doubling model at the GOES angles to the mean 
GOES-derived reflectance of the square area around 
each lidar site. The pixels in the same areas (which 
include the strips analyzed earlier) were analyzed with 
the parameterization using the values of Xc. 

A comparison of the results showed that the average 
value of t v for the radiative transfer model using hex- 
agonal crystals was 0.95 compared to an average of 
0.99 for the parameterization. The rms difference was 
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17%. Although the 0.04 mean difference is not statis- 
tically significant, it is expected for several reasons. The 
Takano et al. ( 1989) results include all of the pixels in 
the region, while the results from ( 5 ) exclude all dark 
pixels. Furthermore, the radiative transfer model does 
not include ozone absorption. Exclusion of the dark 
pixels and inclusion of ozone absorption would both 
tend to increase the mean model-derived optical depth 
by a few percent. The rms differences are most likely 
due to the use of a Lambertian surface albedo in the 
adding-doubling model and the use of a fixed value of 
Xc for each hour in (5). The mean value of £ for the 
model results was 2.55 compared to 2.56 for the pa- 
rameterization. From these comparisons it is evident 
that (5) is quite adequate for parameterizing the de- 
tailed radiative transfer model for the range of optical 
depths seen in the case study. In addition, the good 
agreement between the average derived values of £' over 
land and water at the same hour (i.e., 1600 UTC) in 
Table 4 indicates that ( 5 ) is useful over a reasonable 
range of surface albedos. 

The cirrus clouds observed during the IFO contained 
mixtures of hexagonal columns, plates, and aggregates 
of various sizes and shapes ( Kajikawa and Heymsfield 
1989). Even water droplets were found in some of 
the warmer cirrus clouds (Sassen et al. 1990). In this 
study, the visible radiances have been treated as if they 
were reflected from a cloud composed entirely of one 
size of hexagonal columns. Although exact evaluation 
of the uncertainty in r v due to the use of the single 
phase function is not possible, some inferences may be 
drawn from the results to make a qualitative assessment 
of the hexagonal column, cloud-albedo model. 

One indication of the reasonability of the model used 
here is that the mean value of £ is close to the limit of 
2.0 expected for large particles, especially for the 
warmer clouds. An examination of the results of Tak- 
ano et al. ( 1989) reveal that the mean value of t v de- 
rived with the equivalent ice-sphere model was 30% 
greater than that for the hexagonal column model. The 
resulting mean value of £ for the ice spheres is 2.9. 
These results indicate that the ice-crystal model is more 
appropriate than the spherical model for interpreting 
the cirrus reflectance data. 

The adding-doubling model of Takano and Liou 
( 1989b) was used to calculate the cloud reflectances 
for the angles corresponding to GOES and AVHRR 
in Table 5 using the values of r' v and the hexagonal 
column phase function. The resulting reflectances at 
the GOES angles were ratioed to the corresponding 
reflectances at the AVHRR angles. For the three cases 
in which the GOES and AVHRR angles differed sig- 
nificantly (30 October and 2 November), the average 
ratio for the adding-doubling results was 1.6 compared 
to 1 .3, the mean ratio of the GOES and AVHRR values 
of Xc. About half of the difference in the mean ratios 
may be due to the uncertainties in the values of Xc for 
each individual case. Thus, the actual difference in the 


ratios may be as small as 0. 1 5. This limited result sug- 
gests that the actual scattering by the cirrus clouds may 
be less anisotropic than expected for the hexagonal 
crystal model. Although discrepancies between the mi- 
crophysical characteristics of the actual and modeled 
clouds are likely to cause differences between the ob- 
served and modeled radiation fields, some of the ten- 
dency for more diffuse reflectance from the cirrus may 
be due to shadowing effects as discussed below. A more 
complete analysis using muc! more data is required 
to better define typical cirrus reflectance anisotropy and 
the single scattering phase function that accounts for 
it. Given the current knowledge of cirrus cloud radia- 
tive properties, however, the cloud albedo model used 
in (5) appears to be adequate for computing r' v . 

d. Dark pixels 

It is possible that the “dark” pixels noted earlier may 
be caused by the extreme reflectance anisotropy ex- 
pected for very thin clouds. For small values of r v and 
0 ~ 120°, there should be very little reflection from 
the cloud in the direction of the satellite as compared 
with its albedo to upwelling diffuse radiation from the 
surface. Thus, the pixel may appear darker than a clear 
scene for the same viewing and illumination conditions. 

Examples of the variations of p as a function of a c 
computed with (5) are shown in Fig. 22a for WAU at 
1400 and 1800 UTC. The corresponding values of p s 
are 0.13 and 0.16, respectively, with a ^ = 0.1 1. The 
relationship between cloud albedo and optical depth 
depends on the solar zenith angle resulting in a diver- 
gence of the curves at larger values of a c . Additional 
calculations were performed for po = 0.4; X« = 1; and 
p, = 0.04, 0.08, 0.16, 0.32, and 0.64; with a ^ = 0.92 p,. 
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Fig. 22a. Theoretical top-of-the-atmosphere visible reflectance over 
WAU for 1400 UTC (p, = 0.13) and for 1800 UTC (p, = 0.16). 
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FIg. 22b. Theoretical top-of-the-atmosphere visible albedos for 
various surface albedos. 


The results are shown in Fig. 22b. For relatively dark 
surfaces, the measured reflectance increases monoton- 
ically with a c . Over bright surfaces, the reflectance ac- 
tually decreases first, then increases at greater cloud 
albedos. A dark surface contributes little to the upward 
flux, while the opposite is true for bright surfaces. The 
increased reflectance of a cloud to diffuse radiation 
tends to trap some of the reflected radiation from the 
surface. For thin clouds over bright surfaces, this effect 
can result in a reflectance that is lower than that ob- 
served for the clear-sky case. Additional calculations 
for other solar zenith angles indicate that the insertion 
of the thin cirrus cloud above the bright surface de- 
creases the albedo at lower sun angles, making the dis- 
crimination of clear and cloudy skies more difficult 
over deserts and other bright scenes. It is possible that 
this effect may also be important at certain angles and 
values of X* over dark surfaces. 

Although such a mechanism may cause some pixel 
darkening, it is probably not the primary cause. Cloud 
shading is more likely to be the predominant cause of 
dark pixels. For example, at altitudes of 10 km, a 2 
km thick cloud can cast a shadow that has its leading 
edge removed at least 8 km horizontally on the surface 
from the position of the cloud’s leading edge for do 
» 45°. Thus, the cloud and the portion of the surface 
affected by direct transmission through the cloud can 
easily be offset by one or more GOES pixels. The pro- 
jected line of sight from the satellite through a cloud 
will not line up with the cloud’s shadow unless 0 
«s 180°. In most instances, therefore, the direct trans- 
mission from the surface through the viewed cloud will 
come from the shadow of a different cloud. The diffuse 
reflectance will result from the combined effects of the 
various clouds in the vicinity. For a homogeneous 


cloud field, an optically thick viewed cloud, or 0 
= 180°, such effects are negligible. Figures 1, 2, and 
17 provide ample evidence that the clouds observed 
during the case study are neither optically thick nor 
homogeneous. In an inhomogeneous cloud field, it is 
possible to observe a surface darkened by the shadow 
of a thick cloud through a thin cloud that produces 
little scattering of its own in the direction of the satellite. 
Thus, to the viewer in space a pixel appears to be a 
cloud by virtue of its cold temperature, but is darker 
than expected for a clear scene. 

This effect may be examined quantitatively by con- 
sidering (5). Use of that model implicitly assumes a 
homogeneous cloud field. To consider an inhomoge- 
neous cloud field, let r, replace r v in (6), where r, is 
the VIS optical depth of the cloud casting a shadow 
onto the surface in the line of sight of the observed 
cloud that has optical depth r v . Since the diffuse ra- 
diation emerging from the bottom of the clouds and 
reflecting from the surface is coming from all of the 
clouds in the field, let etc in the third term of (5) be 
replaced by the average albedo of the observed and the 
shading cloud. The results of solving this inhomoge- 
neous version of ( 5 ) for two different observed clouds 
are displayed in Fig. 23. The circles represent r, = r v . 
Viewing angles were fixed at 0 = 52° and X« varied 
with Bo, which changed with time as in the lower halves 
of Tables 3 and 6. Gear-sky reflectance is denoted 
with the dashed horizontal line. 

The greatest effects of cloud inhomogeneities are seen 
for t» = 0.25. Both shading clouds cause dark pixels 
for Bo < 81°. Some dark pixels occur for r, = t v for 6 0 
< 60°. Calculations using smaller no's produced no 
dark pixels. The thickness of the shading cloud becomes 
more important as no increases. The impact of shading 
on the observed cloud having r v = 1 is less pronounced 
with little likelihood of dark pixels for the range of no 


L J 

o 

z 

< 

h* 

o 


a 

2 

tx 

UJ 

CO 

CD 

O 


0.4 

~ ■ 







$ m 


t_ - SHADING 

CLOUD 




■ 





0.3 

= 4.0 

l 








• 

.T 

= 0.25 





♦ 

{ 



0.2 

-° o 

a 

O 

.....o. 

m 

♦ 

r 

1.0 


® 0 




Or. = 

0.25 

0.1 

o 

o 

o 

□ 

O-T. = 

1.0 





<0 

O t . = 

4.0 

on 

1 

i 

i 

1 

1 


0.1 0.2 

0.3 

0.4 

0.5 

0.6 



COSINE 

SOLAR 

ZENITH 

ANGLE 



FIG. 23. Effect of cloud shading on observed reflectance. Optical 
depth of observed cloud is r v . Solid symbols, r v = 1. Open symbols, 
r v = 0.25. Dashed line is clear-sky reflectance. 
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considered. Shading by the thin cloud (r, = 0.25) ac- 
tually increases the observed reflectance for all cases 
using r v = 1 . As t v increases to larger values, the impact 
of r, will become negligible. This lack of a shadowing 
effect for thicker clouds suggests that the observed re- 
flectance may be biased toward a value that is lower 
than expected when the scene contains both optically 
thick (r v > 2 ) and thin clouds, even if no dark pixels 
occur. 

Shadowing effects can be observed in the data shown 
in Fig. 18. Dark pixels corresponding to the thinnest 
part of the cloud are seen before 2027 UTC. The next 
few pixels to the right are detectable but r v is much 
less than r e . The VIS optical depth for the pixels cor- 
responding 2040 UTC is more than three times the 
value of r e , suggesting enhanced reflectance. In this 
figure, the sun would be located approximately to the 
right of the figure at 0 O = 68°. With a cloud top near 
1 1 km, the thickest part of the cloud would cast a 
shadow approximately 28 km to the left or to the po- 
sition corresponding to 2025 UTC. The view from the 
satellite is at a 35° angle, not perpendicular, to the 
solar plane. Thus, the surface in the line of sight of the 
cloud at 2025 UTC would correspond to the 2020 UTC 
position. That surface would be shaded by a cloud with 
r e as 0.5 orr v « 1 .0 assuming that Fig. 1 8 is an ac- 
curate cross section of the cloud. A dark pixel would 
be expected at 2020 UTC based on the results in Fig. 
23. The first visibly detectable cloud coincides with r e 
~ 0.3 at 2027 UTC. Its value of r v is much less than 
the expected value of 0.6. 

While the shading can explain the darker pixels, it 
does not account for the excessively bright ones at 2040 
UTC. For these bright pixels, it is instructive to refer 
to the cloud structure in Fig. 1 . The vertical thickness 
of the cloud changes rapidly so that the pathlength of 
the incident solar radiation through the cloud is not 
necessarily hsecd 0 . The pathlength may be substantially 
increased because of the cloud structure. For example, 
the sun angle may coincide with the right arm of the 
“V” defining the cloud centered at 2035 UTC in Fig. 
1. This effect would tend to cause a cloud to appear 
brighter than expected for the observed r e . In this case, 
the shadowing and enhancement effects apparently 
cancel to yield a reasonable value of £' for the whole 
scene. 

The occurrence of dark pixels is shown in Fig. 24 as 
a function of < b . Most of the dark pixels are found at 
emittances corresponding to r e < 0.2. Some, however, 
are found for r e = 0.5. Very few dark pixels were ob- 
served near noon when 0 > 160°. Most were found 
in the midmoming and late afternoon when shading 
conditions were favorable. The dark pixels comprise 
only 3% of the data considered here. This percentage 
belies the importance of this effect since shading will 
occur in many instances without producing the easily 
detectable dark pixels. The shading effect will tend to 
reduce the observed reflectance causing an underesti- 



EMITTANCE 

Fig. 24. Percentage of total number of “dark” pixels as a function 
of effective beam emittance, * b . 


mation of the cloud’s thickness and emittance. Cloud 
shadows or their absence will also affect the interpre- 
tation of reflectance anisotropy. Their presence at high 
solar zenith angles will yield lower than normal values 
of r v . 

5. Conclusion 

The cirrus clouds observed during the case study 
days are similar to those observed in previous research. 
By comparison, however, they are not necessarily typ- 
ical of midlatitude cirrus clouds since their bulk prop- 
erties are significantly different from those observed 
during other IFO cirrus days. 

The analysis performed in this study using a com- 
bination of ground-based and satellite instrumentation 
yields some results that are similar to earlier studies 
that utilized other combinations of platforms and in- 
struments. This consistency of results for different ap- 
proaches lends a higher degree of confidence to the 
common findings. In addition, the combined datasets 
used here have provided some valuable new insights 
into the problem of determining cirrus cloud properties 
using VIS and IR radiance pairs from satellites. 

The determination of cloud-top emittance or actual 
cloud-top temperature appears to be feasible based on 
the results of this study. Previously, what has been 
termed cloud-top temperature in emittance-adjusted, 
VIS-IR retrievals is closer to cloud-center temperature. 
The emittance ratio and the dependence of cloud 
thickness on cloud-center temperature derived here 
may be utilized in a scheme to estimate cloud thickness. 
This latter parameter will be useful for computing ra- 
diative divergence for observed cirrus clouds. 
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The results indicate that scattering efficiencies are 
greater for colder clouds. This finding holds promise 
for improvements in determining the IR optical depth 
from reflected VIS data. Additional study, however, is 
required to confirm this conclusion. 

From the examination of the reflectance data and 
scattering efficiencies it is concluded that much work 
remains to adequately describe the scattered radiation 
field for real cirrus clouds. Both theoretical and em- 
pirical bidirectional reflectance models should be de- 
veloped for clouds composed of realistic particle shapes. 
As in the ISCCP model, the patterns should be devel- 
oped for various optical depths. Those models will re- 
quire confirmation with the aid of further observational 
data. The effects of cloud particle scattering in IR ra- 
diative transfer should also be examined observation- 
ally. 

Cloud shadows are a problem for the interpretation 
of cloud reflectances, especially for cirrus. The effects 
are not limited to large solar zenith angles. Relatively 
high viewing zenith angles can produce situations that 
result in the viewing of shadows, even for the near- 
zenith sun. It is apparent that analysis of a single pixel 
is most sensitive to shadowing problems. Some of the 
effects may be diminished through averaging over sev- 
eral pixels. Other approaches to deal with shading need 
to be developed. Validations of cirrus scattering cal- 
culations must also consider these effects. The problem 
of shadows is not as important over water because of 
the low surface albedo. Other problems in remote 
sensing of clouds such as partially cloud-filled pixels 
have not been considered here. Future research efforts 
should address these other factors. 
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APPENDIX A 

Low and thick cloud filtering 

The presence of low clouds contaminates the data 
since the analysis is predicated on the presence of only 
one cloud layer (the cirrus layer may contain two or 
more distinct sublayers). The cloud heights are defined 
for the upper-level cloud (generally, z b > 4 km) de- 
tected by the lidar. In some instances, the occurrence 


of low clouds was either noted by the lidar operators 
or evident in the backscatter intensity plots. In other 
cases, the visual observations were not recorded or the 
low clouds were located within the satellite-defined area 
around the site, but beyond the observer’s horizon. 
The datasets containing low clouds must either be 
eliminated or the pixels contaminated by low clouds 
must be removed. Elimination of all datasets contain- 
ing low clouds would substantially reduce the number 
of samples. 

In order to define a simple preprocessing filter, it is 
first assumed that a c oc r v for small optical depths. It 
is also assumed that, to a first approximation, p = T,a c 
+ p t ( 1 — ate) 2 and = p c . Using the definition of £ 
(§ 3e) and (4), it can be shown that 

«b = 1 - exp(— / coc/m), (7) 

where k is the constant of proportionality. Values for 
k are determined in the following manner. For a given 
site and time, q> is computed with ( 1 ) for each cloudy 
pixel having D > A and T < T t - 3 K. The mean 
emittance is then computed for each value of D cor- 
responding to the approximated value of a c . Equation 
(7) is solved for k using each albedo-emittance pair. 
An average value of k is then computed for a specified 
subset of the data. This procedure is applied using T z 
= T c . 

To filter the data, it is assumed that low clouds are 
generally brighter than cirrus, but their emittances de- 
pend on albedo, as in ( 7 ) , in the same manner as cirrus. 
Thus, if a maximum cloud albedo is known for a given 
cirrus emittance, then it may be assumed that any pixels 
brighter than that maximum contain some low clouds. 
The values of those maximum albedos, ow , must be 
estimated here. To determine those maximum values, 
the histograms were first filtered manually to obtain a 
set containing no low clouds. The maximum albedo 
for each temperature ( emittance ) was extracted for ev- 
ery histogram at a given time. Data taken on 22 Oc- 
tober 1986, were used in the analysis. All maxima for 
a given time were compared to determine the greatest 
values for that time. Very little dependence on time 
was found in this comparison. Therefore, the maxima 
from all hours were combined. The resultant values 
were then used to solve (7) for k using a least-squares 
regression fit. The results were averaged to obtain the 

coefficient, k which is used to define the maximum 

cloud albedo for a given emittance. The results of this 
process yield = 2.4. Thus, for a given time, (7) 
may be solved for aw using Aw and the emittance. 
Any pixel with a c > aw + 0.025 is assumed to contain 
low clouds and is eliminated from further processing. 
The 0.025 increment is an estimate of the uncertainty 
in am., resulting from the use of an approximation. 

An example of this filtering is shown in Fig. A1 for 
the data over FMC at 1930 UTC 28 October when low 
clouds were noted in the observer’s log. Crosses rep- 
resent individual pixels. The solid line represents aw 
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FIG. Al. Example of low-cloud filter application over FMC at 
1930 UTC 28 October 1986. All data to right of solid line correspond 
to low clouds. 


computed using k m** • The vertical line, amax = 0.20, 
is used at low reflectances because of larger uncertain- 
ties at those reflectance levels. This value is based on 
a visual inspection of the data that contained low 
clouds. The filtering process eliminates the pixels that 
clearly contain low clouds (i.e., those to the right of 
the line). Some pixels containing both low and high 
clouds, however, may remain in the data. 

All cases having r v > 5.0 were also removed from 
the dataset. This arbitrary limit was imposed to min- 
imize the uncertainties in r e , r v , and subsequent values 
of £. The absolute error in derived optical depth in- 
creases with increasing emittance. The value of T c is 
also more uncertain for thick clouds than thinner 
clouds. Thus, errors in emittance are greater for these 
thick clouds. Their removal does not significantly affect 
the trends in the results. Inclusion of the thick clouds, 
however, greatly increases the noise in the data. No 
case study scenes were affected by the thick-cloud 
elimination. 


APPENDIX B 

Partially cloud-filled pixels 

The use of (7) or (8) implicitly assumes that the 
pixel is completely filled with cloud. While it is likely 
that some of the pixels are partially clear, definitive 
means for determining which are overcast and which 
are partly cloudy are unavailable. Platt ( 1 . , K .^posed 

a bispectral method to discriminate partly cloudy pixels 
from those containing thin clouds. His technique re- 
quires uniform blocks of constant height and thickness 
for detection of partly cloudy pixels. Detection of thin 
clouds allows tHnrVccsses but not ^ 


seen in Figs. 1 and 2, real cirrus clouds are anything 
but uniform and constant in altitude and depth. An 
actual cloud field may also contain broken clouds with 
variable optical depths. The nonuniformity of such real 
clouds blurs the distinction between the linear rela- 
tionship for partly cloudy pixels and the exponential 
dependency for variable thickness overcast pixels. Over 
land the ISCCP algorithm uses a threshold of 6 K colder 
than the clear-sky temperature to determine cloudy 
pixels (Rossow et al. 1988). Because of the low spatial 
variability in. T % over the land areas in this study, it is 
assumed here that pixels colder than T t - 3 K are over- 
cast The data are not used, however, if there is evidence 
of any significant breaks in the ceiling from the lidar 
displays or if the mean IR optical depth is less than 
0.1. Very few breaks were seen over the sites on 28 
October. The cirrus observed on 27 October were 
mostly scattered leading to the selection of only two 
scenes for that day. 
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appendix c 

List of symbols 

particle radius 
Planck function at 1 1.5 pm 
observed and clear-sky visible counts 
cloud thickness 

albedo coefficient, albedo coefficient for 
maximum cloud albedo 
absorption and scattering efficiencies 
cloud-top and -center emittance ratio 
observed and clear-sky 1 1 .5 ^m equivalent 
blackbody temperatures 
cloud-center and -top 1 1.5 pm equivalent 
blackbody temperatures 
clear-sky temperature estimated from sur- 
face data 

air temperatures at altitude z and at the 
surface 
time 

ozone abundance 

cloud-base, cloud-center, and cloud-top 
altitudes 

cloud and diffuse cloud albedos 
clear-sky albedo, clear-sky diffuse albedo 
maxim um albedo for high clouds at a 
given emittance 

extinction, scattering, and absorption cross 
sections 

density of particles in the medium 

vertical emittance, beam emittance 

single-scattering angle 

viewing zenith angle, solar zenith angle 

latitude or wavelength 

COS0, cos 0o 

scattering efficiency ratio 
observed, clear-sky, and cloud visible re- 
flectance 

volume extinction coefficient 
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t c , r v infrared absorption and visible scattering 
optical depths 

r« visible scattering optical depth of shading 
cloud 

T«, T c , T„ ozone, downward cloud, and upward 
cloud transmittances 
<t> longitude 

X„ Xc clear-sky and cloud anisotropic reflectance 
factors 

t relative azimuth angle 
IFO Intensive Field Observations ( 1 9 October- 
2 November 1986) 

IR infrared window ( 1 1 .5 jrm ) 

FMC Ft. McCoy, Wisconsin 
UTC Coordinated Universal Time 
MAD Madison, Wisconsin 
WAU Wausau, Wisconsin 
VIS visible (0.65 /im) 
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